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Five dimensional black hole solutions that describe the QCD crossover transition seen in (2 + 1)- 
flavor lattice QCD calculations at zero and nonzero baryon densities are used to obtain predictions for 
the baryon susceptibility, baryon conductivity, baryon diffusion constant, and thermal conductivity 
of the strongly coupled quark-gluon plasma in the range of temperatures 130 MeV <T < 300 MeV 
and baryon chemical potentials 0 < /is < 400 MeV. Diffusive transport is predicted to be suppressed 
in this region of the QCD phase diagram, which is consistent with the existence of a critical end point 
at larger baryon densities. We also calculate the fourth-order baryon susceptibility at zero baryon 
chemical potential and find quantitative agreement with recent lattice results. The baryon transport 
coefficients computed in this paper can be readily implemented in state-of-the-art hydrodynamic 
codes used to investigate the dense QGP currently produced at RHIC’s low energy beam scan. 
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1. Introduction. After more than a decade of intense ex¬ 
perimental w and theoretical investigation, it is now 
widely believed that the strongly interacting quark-gluon 
plasma (QGP) formed in high energy ultrarelativistic 
heavy ion collisions at RHIC and LHC behaves as a 
nearly perfect relativistic fluid [6]. This surprising fea¬ 
ture of the QGP may be best illustrated by the very 
small value of its shear viscosity to entropy density ratio, 
rj/s ~ 0.2, which has been extracted from comparisons 
between viscous relativistic hydrodynamic modeling and 
experimental data (for a recent review, see [7]). Neither 
ordinary hadronic transport [8] nor weak coupling calcu¬ 
lations [5] have been able to explain such a small value 
for this ratio in QCD, though alternative explanations 
have been proposed [MS] that rely on the appearance 
of a different set of degrees of freedom at the crossover 
transition. 

With the advent of the low energy beam energy scan at 
RHIC, a new region in the QCD phase diagram [HdS] 
corresponding to a hot and baryon rich QGP is being 
explored. The ultimate goal of this effort is to find un¬ 
ambiguous experimental evidence for the existence of a 
critical end point (CEP) [TQIT^. While definite experi¬ 
mental proof of a CEP is still lacking (for a different point 
of view, see [18] ), the observed large degree of collectivity 
displayed by hadrons with low transverse momentum in 
low energy RHIC collisions [H] seems to be compatible 
with the hypothesis that a baryon rich QGP may also 
behave as a nearly perfect fluid [20]. This experimental 
observation, together with the success of viscous hydro- 
dynamic modeling in high energy collisions where baryon 
density effects are small, boosted the interest of the theo¬ 
retical community in applying hydrodynamics to describe 
the dynamical evolution of the hot and baryon dense 
strongly interacting QGP formed in these reactions [2T1 - 
I26j where effects from a conserved baryon current cannot 
be neglected. Given the fundamental role played by hy- 


drodynamical transport coefficients in our understanding 
of the QGP, one expects that the knowledge about the 
novel transport coefficients associated with baryon trans¬ 
port may play a key role in the comparison of hydrody¬ 
namic calculations to the experimental data obtained in 
low energy collisions. 

In this case, the underlying effective model within 
which these transport coefficients are to be computed 
must not only appropriately describe the known equilib¬ 
rium properties of the QGP at nonzero baryon densities 
near the QGD crossover transition m, but also necessar¬ 
ily include from the outset its nearly perfect fluidity. Gur- 
rently, the only consistent framework that is able to deal 
with the required complexity behind the strong coupling 
regime of non-Abelian gauge plasmas and their perfect 
fluid-like features is the holographic gauge/gravity dual¬ 
ity p5II30j . Indeed, a universal prediction of this frame¬ 
work obtained a decade ago is that nearly perfect fluidity 
behavior emerges as a universal property of such strongly 
coupled plasmas and, in fact, r\ls = I/Jtt for a large class 
of gauge theories [31]. Since then, the gauge/gravity du¬ 
ality has been applied to obtain many important qual¬ 
itative insights about the strongly-coupled QGP (for a 
review, see [5^1. 

A few years ago it was realized [33j that holography 
could also be used to shed light on the real time, dy¬ 
namical properties of the strongly coupled QGP near the 
crossover phase transition. The general idea put forward 
in [33] [34] (and also in [35l [36] for a pure glue plasma) 
is that, due to strong coupling effects, many of the prop¬ 
erties of the QGD phase transition may be mimicked 
by black hole solutions of a simple 5-dimensional holo¬ 
graphic model that has the metric, and a scalar field 
(the dilaton), </), as the main dynamical degrees of free¬ 
dom. This type of Einstein-Dilaton bottom-up model not 
only displays the same thermodynamic properties as the 
QGP at zero baryon chemical potential pts = 0 but also 
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has nearly perfect fluidity built-in. Recently, this holo¬ 
graphic setup has been used to investigate, in a quanti¬ 
tative manner, the near-crossover behavior of a large set 
of physical observables (see also the study in the 

Veneziano limit in [SD]). While such models may not be 
directly derived from string theory, a broader view of the 
validity of the standard holographic dictionary may be 
invoked, much like as in the case of condensed matter ap¬ 
plications m, to motivate the calculation of dynamical 
properties of the QGP in a regime that had been other¬ 
wise unaccessible to other nonperturbative approaches, 
such as lattice QCD. 

This is especially true in the case of the baryon rich 
QGP formed in low energy heavy ion collisions. In this 
case, lattice calculations for the temperature (T) and 
baryon chemical potential (/is) dependence of transport 
coefficients may not be available for quite some time and 
holography is a natural framework to perform such cal¬ 
culations. In fact, this was the motivation behind the 
study performed in [52] where a holographic Einstein- 
Maxwell-Dilaton (EMD) model was used to understand 
the energy loss and some thermalization properties of the 
QGP at nonzero baryon density. This model, originally 
proposed in naiio] to study a holographic critical end 
point, has its parameters fully fixed by the equation of 
state of the QGP at fis = 0, as in [33]. However, a strik¬ 
ing finding discussed in detail in |S3] was that, without 
the introduction of new parameters, the model (after the 
necessary revisions to accommodate more recent /is = 0 
lattice data) can also quantitatively describe the available 
lattice data [53] for the equation of state of the QGP near 
the crossover transition with baryon chemical potentials 
up to = 400 MeV. Furthermore, our model prediction 
for the electric conductivity of the QGP at /is =0 [H] 
also gives a good description of recent lattice data [54] 
for this observable (see Fig. 6 in [55] for a direct compar¬ 
ison). These findings give a strong indication that this 
holographic model may be used to make quantitative pre¬ 
dictions for the out-of-equilibrium properties of the QGP 
at moderate values of the baryon chemical potential. 

Here we employ the EMD holographic model [33 1331 
152] to make the first realistic predictions for the T and hb 
dependence of the baryon susceptibility, baryon conduc¬ 
tivity, baryon diffusion constant, and thermal conduc¬ 
tivity of the strongly interacting QGP at the crossover 
transition. The baryon susceptibility is found to increase 
with fj,B while the baryon conductivity is somewhat ro¬ 
bust to the presence of a nonzero ^b, though it displays a 
nontrivial temperature dependence. On the other hand, 
the thermal conductivity increases with fiB- Moreover, 
we predict that the baryon diffusion constant is reduced 
with increasing baryon densities, which indicates that dif¬ 
fusive transport in the strongly-coupled QGP is overall 
suppressed in a baryon rich environment. We also cal¬ 
culate the fourth-order baryon susceptibility at = 0 
and And good quantitative agreement with recent lattice 


data |56j near the crossover transition. 

2. Holographic model. The EMD action reads 
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where the gravitational constant the dilaton potential 
V{(j)), and the Maxwell-Dilaton coupling f{(j)) shall be 
dynamically fixed by lattice data at ^b = 0. 

We look for charged black holes of the form 

ds^ = \—h{r)dt'^ + dx^] + y-yy, (2) 

h(r) 


(j) = (/)(r), and A = A^dx^ = $(r)(it, where the radial 
position of the black hole horizon is given by the largest 
root of the equation h{rH) = 0 and the boundary of 
the asymptotically AdS^ space is at r —>■ oo. In [33] it 
was discussed how to obtain numerical solutions for the 
functions A{r), h{r), and $(r) once V{(j)) and /(^) 
are specified. The equations are solved with a pair of 
initial conditions corresponding to the values of and 
the derivative of the gauge field at the horizon. Different 
choices for this pair of initial conditions translate into 
different thermodynamical states of the plasma. The far 
from the horizon, ultraviolet asymptotics are given by 
[33l3n]- ~ Q^(^)) ^(^) ~ ^o'^) ^{t') ~ and 

$(r) « -b with a(r) = -b 

and z/ = 4 — A, where A is the scaling dimension of the 
operator dual to (j). The ultraviolet coefficients 

and (pA enter in the expressions for the gauge field 
theory observables, such as the temperature, the baryon 
chemical potential, the entropy density, and the baryon 
charge density, respectively: 


T = 


A 




Pb 
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where A « 831 MeV is the energy scale conversion factor 
needed to express in physical units the quantities com¬ 
puted using the black holes (this scale is obtained via 
matching the dip of the speed of sound in the model with 
that found on the lattice [33 [33] ). The dilaton potential 
and the gravitational constant are determined by a com¬ 
parison to lattice data for the (2-bl)-flavor QCD equation 
of state at pb =0 [33] while the Maxwell-Dilaton gauge 
coupling is fixed by the lattice data for the baryon sus¬ 
ceptibility at pb = 0 [37]. The results [33] are = 12.5 
and 


Vifi) = -12cosh(0.606^) -b 

sech(1.2(;i-0.69) 2e-i°°'^ 


/('(') = 


3sech(0.69) 


0.0034(/)®, 

(4) 
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Using Q we numerically generated ^ 10® black hole 
solutions modeling a hot and dense plasma for several 
values of T and ^b- It was shown in |52j that this 
holographic model is in very good agreement with the 
lattice data [S31 for 130 MeV < T < 400 MeV and 
0 < Mb < 400 MeV. At larger temperatures the holo¬ 
graphic model is not applicable for QGP phenomenology 
(since it is not asymptotically free) and weak coupling 
calculations provide a better description of lattice data 

[SS]. 

3. Model predictions. The n-th order baryon susceptibil¬ 
ity is given by Xn = d'^p/dp'% = 9”“V/'9Mb~^- 
we show our prediction for Xa at mb = 0 compared 
with lattice data from |56| (in the same figure we show 
our fit to the lattice data for at mb = 0, which 

was used to determine /(0) in Q). One can see that the 
model can correctly describe the departure of Xa from 
X 2 jT'^ for T’ > 150 MeV, which has been interpreted as 
a signature of deconhnement [59) : at lower temperatures 
in the hadronic gas phase the holographic result for Xa 
is below the lattice data, which illustrates the domain of 
validity of the holographic approach. 
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FIG. 1: (Golor online) Second and fourth-order baryon 
susceptibilities at mb = 0 in our holographic model 
compared with lattice data from |56j . 
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Our predictions for xf different values of mb are 
shown in the top left plot in Fig. One can see that the 
baryon susceptibility increases with increasing mb- This 
behavior can be later verified by lattice QCD calcula¬ 
tions. We also point out that predictions for the isother¬ 
mal compressibility, = p {dp/dp)j., can be directly 
obtained using our results for xf and p via bt = xf / 

m- 

We employ linear response theory to compute the 
baryon conductivity in a holographic setting [ST]. One 
considers linearized perturbations around the back¬ 
ground Q and calculates retarded 2-point functions 
of gauge and diffeomorphism invariant combinations of 
these perturbations, which are then associated with 
transport coefficients via Kubo formulas. As discussed 
in [30], due to spatial isotropy each spatial component 


of the Maxwell perturbation, a^, i = x,y,z, satisfies the 
same decoupled equation of motion and we focus on the 
i-direction (the same result will also hold for the y and 
^-directions), a = a^. Taking a plane wave Ansatz for 
the perturbation with frequency w and zero spatial mo¬ 
mentum, which is sufficient for calculating the DG con¬ 
ductivity (i.e, the w —>■ 0 limit of the spatially uniform 
conductivity tensor), one obtains the equation of motion 

m 


a"(r) + ( 2A'ir) + ) a'(^) 


o-2A(r) 


h{r) \h{r) 


- /(</>)®'W^^ a(r) = 0. (5) 


The Kubo formula for the DG conductivity in phys¬ 
ical units is (excluding the delta function present in 
translationally-invariant systems at finite density [62) 1 


(Tb = - 


A 


lim 

cj —>-0 


[a*{r, uj)a'{r, w)] 

LU 

( 6 ) 


where the expression [a*a'] corresponds to 

a radially-conserved flux such that (§ may be evalu¬ 
ated at any value of the radial coordinate, r. In §, 
following the prescription proposed in [63] to calculate 
retarded propagators in holography, we consider on-shell 
configurations of the vector perturbation corresponding 
to in-falling modes at the black hole horizon which are 
normalized to unity at the boundary. Then, one may 
solve ([^ for a large grid of initial conditions, plug the so¬ 
lutions with the above specihed boundary conditions in 
the Kubo formula and obtain the baryon DG conduc¬ 
tivity as a function of temperature and baryon chemical 
potential. 

The corresponding predictions for this transport coef¬ 
ficient are displayed in the top right plot in Fig. We 
find that the dense QGP at the phase transition is not a 
good baryon conductor. For instance, in the region of the 
(T, Mb)— plane considered here, the ratio T^ctb/?? < 0.05 
in our calculations, which is roughly an order of magni¬ 
tude lower than recent kinetic theory results [64]. One 
can see that ub/T varies signihcantly with T in this 
regime and that a simple conformal approximation (the 
straight line in this plot) would completely miss the im¬ 
portant effects of the crossover transition. One can also 
see that above T ^ 190 MeV the baryon conductivity is 
slightly reduced as one increases the chemical potential 
with the opposite behavior being observed below that 
temperature (the same qualitative behavior was found in 

m)- 

More importantly, the overall dependence of obIT 
with Mb is relatively small even though this class of 
models displays a GEP at larger densities. This is in 
agreement with the fact that holographic models gener¬ 
ally correspond to the type B dynamic universality class 









4 






FIG. 2: (Color online) Baryon susceptibility (top left), baryon DC conductivity (top right), baryon diffusion constant 
(bottom left), and thermal conductivity (bottom right) as functions of the temperature for different values of the baryon 
chemical potential. The corresponding conformal values attained at large T (and = 0) are also shown. 


[55] as pointed out in [55] and, thus, this transport co¬ 
efficient should remain finite as one approaches a CEP. 
Even though QCD has been argued to belong to the type 
H dynamic universality class m, for the values of 
considered in this work (that are not yet in the critical 
region) there may not be much difference between H and 
B classes [55| when it comes to the fluid properties of 
the system. In fact, the expected dynamical critical ex¬ 
ponent for the shear viscosity in QCD ^ 0.054 m 
is small, which will not spoil the perfect fluidity of the 
QGP in low energy collisions. Thus, holographic models 
in which 77 /s does not show singular behavior |66| may be 
applicable to study the perfect fluid transport properties 
of the QGP at not so large baryon densities (outside the 
critical region). In fact, we note that X 2 diverges at the 
CEP while the baryon diffusion constant (to be discussed 
below) vanishes for both H and B classes. 

As shown in [68] . the baryon diffusion may be calcu¬ 
lated using the Nernst-Einstein’s relation, Db = crslxi ■ 
Our results for this transport coefficient are shown in the 
bottom left plot in Fig. We note that baryon diffusion 
already at /is = 0 is predicted to vary significantly across 
the crossover transition approaching its conformal limit 
[61] only at large temperatures. Furthermore, this plot 
shows that there should be a significant suppression of 
baryon diffusion in a baryon rich plasma, even before one 
reaches a (putative) CEP. This precocious suppression is 


a consequence of the robustness of ctb to the presence of 
a nonzero and the enhancement of X 2 ■ 

The bottom right panel in Pig. shows our results 
for the thermal conductivity at the crossover transition 
for several values of fiB- Following m. this may be 
calculated as At = (ctb/T) [(e +p)/p] . One can see that 
this transport coefficient increases with pB and that its 
conformal limit [69] is only reached at temperatures T ^ 
300 MeV. 

Finally, we point out that in hydrodynamic applica¬ 
tions a convenient expression for the dissipative current 
is J = —kbV(p,b/T) with kb = Tgb [70] and, thus, the 
results for the baryon conductivity presented here can 
be readily used to investigate baryon transport effects on 
the QGP in realistic hydrodynamic calculations (tabu¬ 
lated data for gb{T,pb) follow attached as an ancillary 
file). 

4. Conclusions. In this work, the baryon transport 
properties of a baryon rich QGP near the crossover 
phase transition were computed holographically using 5- 
dimensional black holes which quantitatively mimic the 
thermodynamics of (2-1-1)-flavor lattice QCD data at zero 
and nonzero baryon chemical potential [52]. We made 
predictions for the baryon susceptibility, the baryon DC 
conductivity, the baryon diffusion constant, and the ther¬ 
mal conductivity of the strongly coupled quark-gluon 
plasma in the range of temperatures 130 MeV < T < 
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300 MeV and baryon chemical potentials 0 < < 

400 MeV. Baryon diffusive transport is predicted to be 
precociously suppressed in this region of the QCD phase 
diagram even though a critical end point could only ap¬ 
pear in this model at much larger baryon densities. We 
also calculated the fourth-order baryon susceptibility at 
zero chemical potential and found quantitative agreement 
with recent lattice results around the crossover region. 
The transport coefficients computed in this paper can 
be readily implemented in state-of-the-art hydrodynamic 
models used to investigate the dense QGP currently pro¬ 
duced at RHIC’s low energy beam scan or at the CBM 
experiment at the future FAIR facility at GSI. 
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